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Abstract 

We estimated physical conditions of molecular gas in the central molecular zone (CMZ) of the Galaxy, 
using our CO J=3-2 data obtained with the Atacama Submillimeter Telescope Experiment (ASTE) in 
conjunction with J=l-0 12 CO and 13 CO data previously observed with the NRO 45m telescope. The large 
velocity gradient (LVG) approximation was employed. Distributions of gas density, kinetic temperature, 
and CO column density are derived as functions of position and velocity for the entire coverage of the CO 
J=3-2 data. We fairly determined physical conditions for 69% of data points in the CMZ with > la CO 
detections. Kinetic temperature was found to be roughly uniform in the CMZ, while gas density is higher 
in the 120-pc star forming ring than in the outer dust lanes. Physical conditions of high J=3-2/J=l-0 
features are also discussed. 
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1. Introduction 

Molecular gas concentrates in the inner r < 200 pc of 
the Galaxy, the Central Molecular Zone (CMZ). The av- 
erage density and kinetic temperature are much higher in 
the CMZ than anywhere else in the Galaxy (e.g., Morris & 
Serabyn 1996). The distribution and kinematics of molec- 
ular gas there is highly complex. Its small-scale structure 
is characterized by arcs/shells, filaments, and a popula- 
tion of high- velocity compact clouds (HVCCs), suggesting 
a number of supernova explosions may be responsible for 
the boisterous gas kinematics and peculiar physical con- 
ditions there. 

We are performing a large-scale, high-resolution 
CO J=3-2 survey of the CMZ with the Atacama 
Submillimeter Telescope Experiment (ASTE) to pick up 
shocked gas in the CMZ completely (Oka et al. 2006). 
As expected, we have detected many high J=3-2/J=l- 
features in the Al x Ab = 2° x 0.5° area of the CMZ. 
To understand the nature of these high ratio features, it 
is crucial to know physical conditions in them. This pa- 
per presents our attempts to estimate physical conditions 
of molecular gas using our CO J=3-2 data in conjunction 
with J=l-0 12 CO and 13 CO data previously obtained with 
the NRO 45m telescope (Oka et al. 1998). The distance 
to the Galactic Center is assumed to be 8.5 kpc. 

2. The Data 

We used three data sets of CO rotational lines. The first 
data set is the CO J=3-2 line obtained recently with the 
ASTE (Oka et al. 2006). The others are the CO J=l-0 
and 13 CO J =1-0 data obtained with the Nobeyama Radio 
Observatory 45 m telescope (Oka et al. 1998). These data 



sets cover a common region —0.60° < I < +1.50°, — 0.23°< 
b < +0.23°. The velocity coverage is from V LS r = -200 
km s _1 to 200 km s _1 . These data were regrided onto a 
34 // x34 // x2 km s _1 grid, and smoothed to a l / xl / x2 km 
s _1 (FWHM) with a Gaussian filter function. Resultant 
1 a noise levels in the T m b scale are summarized in Table 
1. Total number of data points in this grid is 2284800. 

Table 1. Data sets 



Line Telescope Noise level Reference 

CO J=3-2 ASTE 0.43 K Oka et al. (2006) 

CO J=l-0 NRO 45 m 0.32 K Oka et al. (1998) 

13 CO J =1-0 NRO 45 m 0.28 K Oka et al. (1998) 



3. Analysis 

The large velocity gradient (LVG) model (Goldreich & 
Kwan 1974; Scoville & Solomon 1974) is employed to cal- 
culate intensities of line emission from a small volume 
within a molecular cloud. We calculated CO line intensi- 
ties as functions of CO column density per unit velocity 
width Nco/dV, kinetic temperature Tk, and number den- 
sity of molecular hydrogen n(H.2) (hereafter we write it n). 
We developed a numerical method to derive the physical 
conditions (Nco/dV, Tk, n) from given line intensities of 
CO J=l-0, 13 CO J=l-0, and CO J=3-2 emission (Ti_ , 
T13, Ts-2)- This calculation is a least square fit of the 
LVG model to the observed intensities. Here we describe 
our model, the least square fit, and its application to the 
data. 
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3.1. LVG model 

Our LVG model uses the escape probability of a spher- 
ical cloud. The rotational levels of CO up to J=40 are 
considered. The model uses the collisional rate coefficients 
for de-excitation calculated by Flower (2001) and extrap- 
olated by Schoier et al. (2005). Collision partners are H2 
molecules, and the ortho/para ratio was assumed to be 
3. The 12 CO/ 13 CO abundance ratio is taken to be 24 
(Langer & Penzias 1990). We did not make any explicit 
assumption on velocity gradient, CO/H2 abundance, or 
depth of emitting region along the line of sight. 

The equations of statistical equilibrium are solved it- 
eratively. Physical conditions on a grid which covers 
a domain T>, 10 14 cm~ 2 (km s -1 ) -1 < N C o/dV < 10 19 
cm- 2 (km s -1 )- 1 , 5 K < T k < 120 K, and 10 1 cm -3 < n < 
10 7 cm -3 . This parameter range can cover whole observed 
line intensities above the detection limit. The grid spac- 
ings are 0.1 dex for Nco/dV and n, 1 K for Tk. For 
off-grid points, we analytically interpolate the line inten- 
sities on the grid, using a tricubic interpolation (Lekien 
& Marsden 2005 for general discussions), to reduce the 
computational cost. 

3.2. Fitting Procedure 

The problem we have to solve is to determine 
x={N C o/dV, T k , n) for given j/=(Ti_ , ^3^3-2), whose 
elements are observed line intensities of CO J =1-0, 13 CO 
J =1-0, and CO J =3-2 respectively. For this purpose, we 
employed the least-square fit method. The best fit param- 
eter x gives a local minimum of the chi-square function 

X \x) = Y J {yi-T;{x)¥/al (1) 

where T*(x) is the line intensity predicted by the model, 
and G{ is 1 a noise level of y\. Note that the function 
is known as chi-square distribution with 3 degrees of free- 
dom. The chi-square function become zero when the given 
line intensities y can be reproduced by the model. We can 
assess the quality of fit x by the value of x 2 . 

Line intensities converge to ~ Tk in the optically thick 
thermalized region. In this region, it is difficult to con- 
strain n and Nco/dV, since line intensities stay nearly 
constant over wide range of parameters. A simple min- 
imization of x 2 often returns unphysically high n and 
Nco/dV in the thermalized region. To suppress this di- 
vergence, we introduced a constraint function J cons (a;) and 
searched the minimum of x 2 ( x ) +/consW- The constraint 
function is written as 

{Ccons ( x ^ D) 

(*6A»<0.5»£ t a ) (2) 
\ ' crit / 

where C CO ns is an arbitrary large constant to restrict the 
solution to domain D, n 3 ^" 2 is the critical density of the 
CO J=3-2 transition (10 3,8 cm -3 ), and a is a positive 
constant which determines the strength of modification 
to suppress unrealistic high density. The form of / cons (x) 
on D is chosen so that it is in n < 0.5n 3 ^" 2 , quadratic 
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in n > 0.5n 3 ^" 2 , differentiate on D, and a is taken so 
that /cons = 1 at n = n 3 T 2 . Figure 1 shows an example of 

X 2 0) + /cons (^) • 

To solve the minimization problem, we tested the steep- 
est descendent method and the quasi-Newton method 
BFGS algorithm. Although they gave the same result 
in most cases, the latter is less stable and sometimes 
failed to get a global minimum. Thus we employed the 
former method in our calculations. We used a fixed ini- 
tial value 7V C o/dV=10 16 cm -2 (km s" 1 ) -1 , T k =25 K, and 
n= 10 3 cm -3 , for all parameter range. A local minimum 
of x 2 ( x ) + icons (^) was obtained in any case. We used 
X 2 to evaluate the propriety of solutions. If x 2 — H-35, 
which corresponds to 99 % confidence of the chi-square 
distribution with 3 degrees of freedom, we can say that 
the parameters are well-determined. 

3.3. Existence Examination of Solution 

Figure 2 shows results of the least-square fit for T3_2 = 
10 K and various Ti_o and T13. Parameters are well- 
determined in a region shown in gray. Vertexes on the 
boundary of the well-determined region are named A, B, 
C, and D. Each part of the boundary corresponds to a lim- 
iting case. Segment AB corresponds to the optically-thick 
limit where Ti_o = Ti3 C0 . Segments BC and CD corre- 
spond to the local thermodynamical equilibrium (LTE) 
case for optically thick and thin limits respectively. This 
boundary comes from a condition that excitation tempera- 
ture of J =1-0 is higher than that of J=3-2. Segments DE 
and EA correspond to the subthermal limit for optically 
thin and thick respectively. Kinetic temperature increases 
as the excitation temperature of J=3-2 becomes subther- 
mal when lines are optically thin, and density decreases 
when optically thick. Out of the well-determined region, 
there is a line of discontinuity in parameters, indicated by 
a bold line in Figure 2. 

3.4- Application to the Data 

We applied the fitting method for each pixel of the 
data described in §2. We assume that beam filling factors 
(r/i_o, r/13, 773-2) are 1 for all lines. The smaller beam fill- 
ing factors, 77 < 1, generally give a higher temperature for 
most data points. In case of 7713 < 771-0 = 773-2, which of- 
ten occurs in giant molecular clouds, kinetic temperature 
decreases while density increases slightly. 

Another breakdown of one-zone assumption occurs in 
a highly optically thick cloud. In such a cloud, highly 
opaque 12 CO line is severly self-absorbed, tracing the low- 
density, subthermally excited envelope, while less opaque 
13 CO line traces the high-density cloud core. This sit- 
uation gives T13 similar to T12, thereby one-zone LVG 
analyses give low- temperature, high-density solutions. 
Unfortunately, x 2 of sucn improper solutions are not nec- 
essarily high. A discontinuity of physical conditions ap- 
pears at the boundary between the proper and improper 
solution areas in the l-b-V cube. These improper solutions 
can be rejected by the behavior of x 2 around the mini- 
mum; x 2 around improper solutions are more sensitive to 
the temperature than those around proper solutions. We 
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Fig. 1. Function X 2 (x) + /cons(x) for Ti_ =13 K, T 13 =l K, and T 3 _ 2 =10 K, (a) at N C o/dV=10 16 - 3 cm- 2 (km s -1 ) -1 , (b) at 
n=10 3 ' 55 cm -3 . Contours are 11.35 and 50. The cross indicates the minimum of the function. 



employed the rejection criterion, d 2 x 2 /dT^ > 10. 

The LTE assumption may hold well in such thermalized 
optically thick cloud cores. In this case, T13 can be related 
to TVco by, 

/T 13 dy=(/ 12/1 3 + l)-Iexp(-^-) 

T13 = {/.(T ex ) - J„(T bg )}(l - e~ T ) 

where A2/13 is the abundance ratio 12 CO/ 13 CO, Z is the 
partition function, v is the frequency, T ex is the excitation 
temperature, A\§ is the Einstein A-coefflcient, I V {T) is the 
Planck function, Tb g is background radiation temperature, 
and r is the line optical depth. Reducing eq.(3), we obtain 

7V C o/dy = 1.15 x 10 15 (r ex + 2.64exp|^ 

This means that Nco/dV is approximately proportional 
to T13. If T ex varies from 20 K to 50 K, Nco/dV increases 
~ 70 %. The discrepancy within a factor of two may not 
be serious for Nco/dV, considering it ranges over three 
orders of magnitude. 

4. Results and Discussion 

We successfully obtained proper solutions with signifi- 
cant confidence (x 2 < 11.35) for 400841 data points out of 
653550 points. Histograms of physical conditions and x 2 
are shown in Figure 3. The lower limit of Nco/dV is set 
by detection limits of spectral lines. Nco/dV reaches 10 19 
cm -2 (km s -1 ) -1 only for the foreground arms where line 
optical depths are quite large, where the parameters con- 
tain large uncertainties. We had low kinetic temperatures 
(~5 K) in the foreground arms. The highest temperature 



(^120 K) was found in the vicinity of Sgr A. 

4.1. Distribution of physical conditions 

Figure 4 shows longitude- velocity diagrams of physical 
conditions at b=— 0.009° for example. Data in the velocity 
range —60 km s _1 < Vlsr < 20 km s _1 , are severely con- 
taminated by the four spiral arms in the Galactic disk. 
However, since our analyses are inapplicable to opaque 
gas, we exclude this velocity range from the following 
discussion. For the data points in the velocity ranges, 
Vlsr < —60 km s _1 and 20 km s _1 < Vlsr? we have proper 
LVG solutions for 281885 data points out of 410704 data 
points (68.6 % ). Here we describe physical conditions of 
the two major components in the CMZ: the 120-pc star 
forming ring (Sofue 1995) and the outer dust lanes. 

4.1.1. Star- Forming Ring 

A large fraction of dense molecular gas in the CMZ is 
confined in the 120-pc star- forming ring, which contains 
Sgr Bl, B2 and Sgr C HII regions. Our calculations show 
that the star-forming ring typically has kinetic tempera- 
ture of 20 - 35 K and density of IO 3 5 - 4 cm" 3 . There 
are some areas where our LVG model cannot reproduce its 
line intensities, being associated with spots of improper 
solutions. In these areas, the opaque CO J =1-0 line is 
likely self-absorbed by low-density envelope of molecular 
clouds. 

The result is consistent with previous works of NH3 
lines (Hiittemeister et al. 1993) and those of high den- 
sity tracers (e.g. Tsuboi et al. 1999). The low- J CO 
lines arise from the cool, dense component suggested by 
Hiittemeister et al. (1993). Martin et al. (2004) reported 
the results of LVG analysis using [CI], CO J =7-6 and 
J=4-3 lines. They found kinetic temperatures about 50 
K; this is higher than ours, and the discrepancy may be 
due to the difference in transitions used. 

4.1.2. Outer dust lanes 

The structures K, L, and M (Rodriguez-Fernandez et 
al. 2006), which had been recognized as parts of the 200- 
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Fig. 2. Physical conditions determined by the LVG model. The CO J=3-2 line intensity is fixed to 10 K: (a) CO column density 
per unit velocity width, (6) kinetic temperature, (c) number density of molecular hydrogen, and (d) chi-square function. In (a) and 
(c) contours are labeled by logarithm to base 10. A region where chi-square is less than 11.35 is shown in gray. Thick lines indicate 
the discontinuity line of best-fit parameters. 



pc expanding molecular ring (Kaifu et al. 1972), are now 
understood as dust lanes associated with the central bar. 
They have roughly uniform kinetic temperature about 30 
K and low density io 3 0-3 - 5 cm -3 . The density is consis- 
tent with the non-detections of high-density tracers. Hot 
(250 K) and low-density (< 10 2 cm -3 ) component associ- 
ated with the structure M has been detected by the H3" 
absorption line observations toward the Quintuplet cluster 
(Oka et al. 2005). However, we found no strong evidence 
of high temperature in the structure M. Two explanations 
for the disagreement may be possible; (1) hot gas detected 
in H3" absorption has extremely low density which cannot 
excite the J=l level of CO; (2) hot gas is spatially local- 
ized to the Quintuplet cluster. 

4.2. Mass Estimate 

We made a map of CO column density in the CMZ 
by summing up Nco/dV over the velocity ranges from 
Vlsr = —200 km s _1 to —60 km s _1 and from +20 km 
s _1 to +200 km s _1 (Fig. 5). The maximum Nqo is 
3.1 xlO 19 cm -2 found in the direction of the Sgr B2 core. 



This value is consistent with the 13 CO column density 
9xl0 17 cm -2 (Lis & Goldsmith 1989). Assuming that 
[CO]/[H 2 ] is 2.4xl0- 5 , based on [ 13 CO]/[H 2 ] given by Lis 
& Goldsmith (1989), molecular mass included in our data 
is ~ 2 x 10 7 Mq. This amount is smaller than the known 
value of molecular gas mass in the CMZ, 5-10 xlO 7 M 
(Morris & Serabyn 1996). The discrepancy may be ac- 
counted by; (1) incomplete coverage of our data for the 
entire CMZ, (2) inapplicability of our analyses to high 
Nco/dV regions. 

4.3. High CO J =3-2/ J =1-0 Intensity Ratio Gas 

The CO J=3-2/J=l-0 intensity ratio, #3-2/1-0 = 
T3-2/T1-0 may be useful to diagnose the physical con- 
ditions. Generally, high R 3 _ 2 /i-o exceeding 1.5 is a sig- 
nature of highly-excited, optically thin gas with high tem- 
perature and high density. Figure 6 shows the CO J=l- 
distribution of high temperature gas, Tk > 60 K, and 
high density gas n > 10 3,7 cm -3 . The distribution of high 
#3-2/1-0 S as (Fig. 5 of Oka et al. 2006) is strikingly 
similar to that of high temperature gas. This suggests 
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Fig. 3. Frequency distribution of the physical conditions: (a) CO column density per unit velocity dispersion interval, (6) kinetic 
temperature, (c) number density of molecular hydrogen, and (d) chi-square function. Improper solutions with chi-square less than 
11.35 are shown in light gray. Proper solutions with chi-square less than 11.35 are shown in dark gray. 



that i?3_ 2 /i_o is a good measure of kinetic temperature, 
demonstrating validity of extracting shocked gas by high 

^3-2/1-0- 

5. Summary 

We developed a robust scheme to estimate physical con- 
ditions by fitting results of LVG calculations to multi-line 
data sets. We applied the scheme to the CO J=3-2 data 
of the Galactic center CMZ obtained with the ASTE and 
the CO, 13 CO J=l-0 data sets in the same region ob- 
tained with the NRO 45 m telescope. The analyses made 
a success to some degree, properly determining physical 
conditions within 99% confidence level for about 69 % of 
the data points. The physical conditions in the CMZ are 
summarized as follows; 

• the 120-pc star forming ring has higher density than 
the outer dust lanes. 

• the kinetic temperature of the outer dust lanes are 
uniformly low, ~ 30 K, 

• high R 3 _ 2 /i-o (> 1-5) gas in this region is mostly 
high temperature gas. 

We could not determine physical conditions in high- 
density cores with large column density, since we used two 
opaque lines with rather low effective critical densities. 
Breakdowns of the one-zone assumption might also 
worsen the accuracy of our analysis. The use of less 
opaque lines such as C 18 0, as well as high-density tracers 
must improve the accuracy. 
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Fig. 4. Longitude- velocity diagrams at b=— 0.009° of the physical conditions (color); (a) CO column density per unit velocity 
dispersion interval, (b) kinetic temperature, (c) number density of molecular hydrogen, and (d) chi-square function. The diagram 
of chi-square function (gray) and areas of improper solutions (purple) are overlaid in (6) and (c). The structures K, L, and M 
(Rodriguez-Fernandez et al. 2006) are indicated in (a). 
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Fig. 5. A map of CO column density for the velocity ranges from Vlsr = —200 km s _1 to —60 km s _1 and from +20 km s _1 to 
+200 km s _1 . 




Fig. 6. (a) A map of CO J=3-2 emission ( / T 3 - 2 dV) integrated over velocity ranges from V LS r = -200 km s" 1 to -60 km s" 1 
and from +20 km s — 1 to +200 km s _1 for data with kinetic temperatures higher than 60 K. (b) Longitude- velocity diagrams of CO 
J =3— 2 emission integrated over available latitudes (^T3_2)for data with kinetic temperatures higher than 60 K. (c) and (d) Same 
as (a) and (6) respectively but with densities higher than 10 3-7 cm -3 . 
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